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Abstract The transverse shrinkage behavior of early wood 
and late wood tracheids of radiata pine (Pinus radiata D. 
Don) was investigated by the power spectrum analysis. The 
representative cell model shapes before and after shrinkage 
constructed by the analysis revealed that the early wood 
tracheid showed anisotropic shrinkage, although the late 
wood tracheid showed almost isotropic shrinkage. To link 
the macroscopic shrinkage of coniferous wood with the re¬ 
sults obtained by the power spectrum analysis, a two-layer 
model composed of early wood and late wood was adopted, 
and the relation between shrinkage anisotropy and late 
wood fraction was predicted. The results suggested that 
the shrinkage anisotropy depended significantly on the 
mechanical interaction between early and late wood. 

Key words Radiata pine • Cell shrinkage ■ Shrinkage 
anisotropy ■ Power spectrum analysis ■ Two-layered model 


Introduction 

The shrinkage anisotropy of coniferous wood in the trans¬ 
verse direction is believed to be affected by the arrange- 
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ment and transverse shapes of tracheids. It is well known 
that the separated early wood shows anisotropic shrinkage, 
whereas the separated late wood shows comparatively 
isotropic shrinkage. 1 ' 2 Coniferous wood from tropical re¬ 
gions without distinct growth rings also shows anisotropic 
shrinkage. The effect of cellular structure on shrinkage 
anisotropy in wood has been investigated through micro¬ 
scopic observations and theoretical approaches using cell 
models. 3 The causes of anisotropic shrinkage are not always 
clear. An analysis using a representative cell model would 
be useful for investigating the shrinkage anisotropy of 
wood, but it is difficult to establish such a model as there is 
relatively wide variation in the cellular structure of wood. 

Maekawa et al. 4 have applied a two-dimensional power 
spectrum analysis to the determination of the anatomical 
features in the transverse section of wood. They quan¬ 
titatively determined the cell arrangement, which had 
been qualitatively recognized by microscopic observations. 
Furthermore, this analysis was utilized to investigate the 
shrinkage behavior of tracheids. 5 The shrinkage of a 
representative model of coniferous wood agreed well with 
macroscopic shrinkage. 

These results suggested that the power spectrum analysis 
can effectively clarify the relation of cellular structures to 
macroscopic shrinkage of coniferous wood. In this study, 
the transverse shrinkage behavior of early and late wood 
tracheids of radiata pine were investigated using the power 
spectrum analysis. A two-layered model composed of early 
wood and late wood was chosen to link the macroscopic 
shrinkage of coniferous wood with the results obtained by 
the power spectrum analysis, and the relation between 
transverse shrinkage anisotropy and late wood fraction was 
investigated. 


Materials and methods 

Radiata pine (Pinus radiata D. Don) with an air-dried den¬ 
sity of 0.41 g/cnT and an average annual ring width of 15 mm 
was used. Specimens for measurement of shrinkage and 
Young’s modulus were prepared from the same block. As 
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shown in Fig.l, five early wood specimens [A: 5mm (R) X 
5mm (T) X 40mm (L)], five late wood specimens [B: 5mm 
(R) X 5mm (T) X 40mm (L)] and ten specimens containing 
late wood and early wood [C: 30 mm (R) X 30 mm (T) X 
5 mm (L)] were prepared from the wood block (R, radial; T, 
tangential; L, longitudinal). Lines along the radial and tan¬ 
gential directions were drawn on the transverse surface 
of specimen C to measure shrinkage. All specimens were 
soaked in distilled water until saturated. The transverse 
surfaces of specimens A and B in the wet condition were 
smoothed by a sliding microtome. 

Construction of cell model shape 


derived from the binary figures and was shown by a dot 
(diameter 7 pixels) on the monitor. After the dots and back¬ 
ground were converted to a gray-level image of 255 and 0, 
respectively, a fast Fourier transform (FFT) was conducted 
to obtain a power spectrum. 

The periodicity of the cell arrangement was determined 
from the power spectrum as follows 4 : A set of dots is math¬ 
ematically expressed by the following function D(s, t) which 
convolutes the delta functions of cell positions 8{s - x h 
t — y,) and the aperture function of a dot A(s, t). 

D(s, t) =A(s, t)*{<5(r -x u t-y) +... 5(s -x N , t -y w )} 

( 1 ) 


After all surfaces of specimen A, except the one to be 
replicated, had been sealed with a thin polyvinylidene chlo¬ 
ride film to avoid water evaporation, the specimen was a 
pressed against polyethylene film (0.1mm thickness) soft¬ 
ened by heating at about 130°C on a glass slide to obtain a 
replica. 6 The specimen was dried slowly at room tempera¬ 
ture for 5 days, placed in a closed chamber conditioned at 
98% relative humidity (RF1) for 5 days, dried over P 2 O s at 
60°C, and finally dried overnight at 105°C. A replica under 
oven-dried conditions was then obtained. The late wood 
portion of specimen B was separated using a razor. Replicas 
under both conditions were obtained in the same manner. 

Light micrographs were obtained at three positions on 
the replicas prepared under wet and oven-dried conditions. 
Cell matrices of 15 X 15 were set up on the micrographs. 
The edges between walls and lumens in the matrices were 
traced on a clear sheet. The traced figures were then input 
into an image analyzing system (Luzex III) using a televi¬ 
sion camera so the radial direction of figures coincided with 
the axial direction of the monitor, smoothed and converted 
to binary figures. The center of gravity for each cell was 
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5 mm (R) X 5 mm (T) X 40 mm (L) 

Fig. 1. Preparation of specimens for measuring shrinkage. R, radial 
direction; T, tangential direction; L, longitudinal direction 


where x, and y, are coordinates of the zth dot, and * is 
an operator of convolution. By the Fourier integration of 
D(s, f), the power spectrum of the delta functions in polar 
coordinates ( k, (p) is given by 


Pg{k, <p) = 
(**/) 


1 N N 27Efc{(*, -* y )cos (p(yi -y,)sin<p} 

2 hh cos M 


+ N 
( 2 ) 


where k is the frequency, and M is the size of the monitor 
(1024 in this case). This equation represents the Fourier 
transform of all dots projected on a line in the <p direction. 
The periodicity of cell position in the radial direction 
(<p = Till) and the periodicity of radial files in the tangential 
direction (<p = 0) were determined from this equation. 

The direction of cell arrangement can be estimated from 
the power density at each angle for Eq. (2), namely, 


1 N N , 

r.(*HXXC 

L i=\ y=l 1 


2nk< 


(jq —XjJ cos tp 


2 cos 

M 


dk (3) 


The power density of each angle at a low frequency band of 
30-100 pixels was adapted to determine the direction of the 
cell arrangement. The density was obtained by reading and 
integrating the gray levels of the power spectrum image 
using a computer. 

To investigate the cell wall orientation in the transverse 
direction, the cell wall regions between adjacent lumens in 
the figure input were painted by the computer, and a figure 
of the central line of the cell walls was created by thinning 
the thickness to 1 pixel. This figure was converted to a gray- 
level image and subjected to FFT to obtain a power spec¬ 
trum, which was saved on a floppy diskette. The power 
spectrum of the image consisting of the rectangular ele¬ 
ments with length L, width B , and weight per unit area m is 
expressed by 

P r ( k ’ e ) = “=4,f m2 smc 2 {^Lcos(0' - 0)} ^ 

■ sine 2 {nkB sin(0' - • H(m,6')d0'dm 
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where h and m are the number of elements and the weight 
per unit area of the image, respectively, and H (m. O’) is the 
function expressing the relation between m and the longitu¬ 
dinal direction of an element O' , 7 The orientation of ele¬ 
ments can be estimated by the power spectrum density at 
each angle: 

Pr {o) = ['P r {k,e)dk (5) 

The following relation is confirmed by a numerical analysis. 
p r (d + 7t/2) = j H{m,0)dm (6) 

The orientation of the cell wall was estimated by regard¬ 
ing the figure made by thinning the cell walls as consisting of 
rectangular elements. The power density of each angle at a 
higher frequency band of 200-300 pixels was adapted to 
determine the orientation of the cell wall. The density was 
obtained by reading and integrating the gray levels of the 
power spectrum image using a computer. 

The cell model shape was constructed as follows. Cell 
positions were determined from the periodicities and direc¬ 
tions of the cell arrangement. A cell element was placed at 
an angle determined from Eq. (6) in the middle of two cell 
positions. All cell elements were connected to construct the 
model shape, as shown in Fig. 2. 



Fig. 2. Transverse cell model shape constructed by the power spectrum 
analysis. d r , d R , tangential and radial periodicities of cell position deter¬ 
mined by Eq. (2); cp, direction of cell arrangement determined by 
Eq. (3); 6, direction of cell wall orientation determined by Eq. (6). 


Measurement of macroscopic shrinkage and 
late wood fraction 

Specimen C was dried in the manner mentioned above, and 
the radial and tangential shrinkages were calculated from 
the dimensional change of the specimen measured by a 
traveling microscope. The radial width of the darker por¬ 
tion of specimen C was measured by a traveling microscope 
to determine the late wood fraction. 


Measurement of tangential Young’s modulus 

Two wood blocks containing one annual ring [15 mm (R) X 
50 mm (T) X 100 mm (L)] were glued at radial surfaces with 
polyvinyl acetate adhesive. Ten specimens [15 mm (R) X 
100 mm (T) X 1mm (L)] were prepared from the block to 
measure the tangential Young’s modulus. Ten specimens 
[10 mm (R) X 100 mm (T) X 1mm (L)] to measure the 
tangential Young’s modulus of early wood specimen were 
prepared in the same manner. The measurement was car¬ 
ried out at 20°C and 60% RH by a free - free flexural 
vibration method. It has been confirmed that the glue line 
has no influence on the results. 8 


Results and discussion 

Shrinkage behavior of early and late wood tracheids 

Figure 3 shows the representative cell model shapes of early 
and late wood tracheids constructed by the power spectrum 
analysis. The shrinkage determined by the power spectrum 
analysis for the center of gravity of each cell agreed well 
with the macroscopic shrinkage, and the features of shrink¬ 
age deformation could be expressed well by the models 
constructed by the analysis. 5 

The representative model of an early wood tracheid of 
radiata pine had a hexagonal, slightly asymmetrical shape. 
The cell rotated slightly during drying. The radial and tan¬ 
gential shrinkages of the model were 2.5% and 5.5%, re¬ 
spectively. The axial shrinkages of cell elements a, b, and c 
were 4.2%, 3.7%, and 6.2%, respectively. The included 
angles between two radial walls under wet (a) and oven- 
dried (ft) conditions were 128° and 133°, respectively, so 
the angle increased by 5° during drying. The early wood 
tracheid itself showed anisotropic shrinkage with changing 
cell shape, which might cause anisotropic shrinkage in the 
separated early wood. On the other hand, the transverse 
shape of late wood tracheids was almost rectangular. The 
radial and tangential shrinkages of the model were 7.2% 
and 8.0%, respectively. No rotation of cell shape during 
drying was observed. The included angle was slightly in¬ 
creased during drying. The axial shrinkages of the cell ele¬ 
ments d and e were 7.9% and 7.2%, respectively. Late wood 
tracheids showed almost isotropic shrinkage. 

Different shrinkage behavior was observed between 
early and late wood tracheids of radiata pine. Therefore 
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Fig. 3. Cell model shapes under wet (solid line) and oven-dried 
(.dotted line) conditions for early (A) and late (B) wood tracheids of 
radiata pine 


when analyzing the shrinkage anisotropy of coniferous 
wood with a low volume fraction of ray tissue, the mechani¬ 
cal interaction between early and late wood must be taken 
into consideration. 


Shrinkage anisotropy of coniferous wood 


Microscopic shrinkage of early and late wood was deter¬ 
mined by the power spectrum analysis. The radial and tan¬ 
gential shrinkage and the tangential Young’s modulus of 
whole wood specimens and the tangential Young’s modulus 
of early wood specimens were also measured. The results 
are shown in Table 1. Using these values, shrinkage of the 
two-layered model was calculated to investigate the shrink¬ 
age anisotropy of coniferous wood. 

Figure 4 shows a model of coniferous wood in which 
layers a and b represent late and early wood, and the direc¬ 
tions 1 and 2 correspond to the tangential and radial direc¬ 
tions, respectively. Model A under wet condition may 
change to model B after drying. A homogeneous strain may 
occur in direction 1 but not in the edge portions in direction 
2, as shown in model B. The shrinkages [e, (i = 1, 2)] and 
stresses [<S, (i = 1,2)] of the model can be expressed with the 
shrinkages [ef] and stresses [<5f (i = 1, 2; k = a, b)\ of early 
wood and late wood. 

f i = el = el o 1 = (pal + wl (= 0) 

£■2 = 0£?' + W'i c 2 ~ £>i = ( = 0) (7) 

where cp and t// (=1 — <p ) are late wood and early wood 
fractions, respectively. Assuming elastic strains occur in a 
and b layers due to mechanical interaction during drying, 
the shrinkages of a and b layers can be expressed by the 
following equations. 


El 


-a: - - 


a: 


1 




(8-1) 


Table 1. Experimental values for radiata pine 


Tangential shrinkage e, (%) 6.32 

Radial shrinkage % (%) 2.54 

Tangential shrinkage of late wood a" t (%) 8.01 

Tangential shrinkage of early wood a\ (%) 5.46 

Radial shrinkage of late wood ot) (%) 7.23 

Radial shrinkage of early wood a b 2 (%) 2.52 

Tangential Young’s modulus £, (GPa) 0.42 

Tangential Young’s modulus of early wood E\ (GPa) 0.37 



Fig. 4. Coniferous wood model. Layers a and b correspond to late and 
early wood, and directions 1 and 2 to tangential and radial directions, 
respectively. A represents wet conditions and B is after drying 


: —- 2 ^- 0 -“ +al 
El 1 E a 2 2 2 


e b 2=- 


+ J_ 

7 * °2 + 

Ei 


E\ 


a b 2 + a\ 


( 8 - 2 ) 


where Ef (i = 1,2; k = a, b) are the elements of the matrix 
and af (i = 1,2; k = a, b) are the shrinkages of layers. From 
these equations, e, and e, are given by 


= Aa“ + (l - A)ctf 


(9-1) 


e 2 = (pocf + y/af + C 


(9-2) 


£, =<pEl +y/El 


(9-3) 


where A and C are <pE\lE x and (af - af){i//vf 2 A - <pv a n 
(1 - A)}, respectively. The Poisson’s ratios for late and early 
wood portions, vf 2 and vf 2 (the ratio of the strain in the 
direction 1 to that in the direction 2), are required to calcu¬ 
late the shrinkage in the direction 2. Generally, the 
Poisson’s ratio of cellular solids can be determined by geo¬ 
metrical calculation of the cell shape. Regarding the slightly 
asymmetrical shape of the early wood tracheid shown in 
Fig. 3A as a symmetrical one, vf 2 was determined by the 
following expression 9 


(a/b + sin 6) sin 0 


cos" 6 


( 10 ) 


where a/b is the ratio of the axial length of cell element a to 
that of cell element b. The value of vf 2 was calculated to be 
0.74. However, this equation cannot provide the Poisson’s 
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ratio of the late wood portion in which tracheids have lower 
0 values. For convenience, v", was assumed to be 0.5. 

Substituting the experimental values shown in Table 1 in 
Eq. (9-1), a value of 0.25 was obtained for 0, and a value of 
0.57 GPa was estimated for E\. On the other hand, a value 
of 0.06 was obtained by substituting the values for Eq. (9-2). 
In this case, a value of 1.20GPa was estimated for E“. The 
experimental value of 0 for specimen C was 0.17. In this 
case, a value of 0.66 GPa was obtained for E“. There was a 
difference between 0 values calculated from Eq. (9-1) and 
(9-2). This difference was considered to originate in the 
simple analysis using a two-layer model. Elowever, the ex¬ 
perimental values of 0.17 showed an almost middle value 
between the two calculated values of 0.06 and 0.25. There¬ 
fore it was considered that the analysis of transverse shrink¬ 
age behavior by the two-layer model was valid as a first 
approximation. 

The shrinkage anisotropy of coniferous wood was pre¬ 
dicted from Eq. (9-1) and (9-2) with values of 0.17 and 
0.66 GPa for 0 and E\. Figure 5 shows the relation between 
tangential shrinkage £, and late wood fraction 0 . The solid 
line represents the calculated values, which increased with 
increasing 0. Figure 6 shows the relation between radial 
shrinkage and 0. The solid line shows the values calcu¬ 
lated by Eq. (9-2), and the dotted line shows the values 
calculated without the interaction term C. Although both 
values increased with increasing 0 , the dotted line increased 
more rapidly. Figure 7 shows the relation between the 
shrinkage anisotropy (£ 1 / 62 ) and 0. The dotted line de¬ 
creased more rapidly with increasing 0. These results show 
that the influence of the mechanical interaction on £,/% is 
significant. 


Because y/v b u X was about three times larger than 0 v “ 2 
(1 — A), the effect of \j/v b n X on C is remarkable. 


yv b n l = yv[ 2 


<t>E\ + yr£t 


This equation contains early wood fraction (y/) and the 
Poisson’s ratio of early wood (v* 2 ), Young’s modulus of late 



0 


Fig. 6 . Relation between calculated values of radial shrinkage (e,) and 
late wood fraction ( 0 ). Solid and dotted lines show the results calcu¬ 
lated from Eq. (9-2) and the equation without the interaction term C, 
respectively 



Fig. 7. Relation between shrinkage anisotropy (£ 1 / 62 ) and late wood 
Fig. 5. Relation between calculated values of tangential shrinkage (e,) fraction (0). Solid and dotted lines show the results calculated with and 
and late wood fraction (0) without the interaction term C, respectively 
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wood (£“), and Young’s modulus of early wood (£?). E" 
can be regarded as a constant value, because there is almost 
no difference in the transverse shapes of late wood 
tracheids among species. Therefore C depends significantly 
on both v 22 and E\. These results suggest that the shrinkage 
anisotropy of coniferous wood depends remarkably on the 
transverse shape of early wood tracheids. 
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